Activity on the M star of QS Vir. 



Ribeiro, T. 

Departamento de Fisica, Universidade Federal de Santa Catarina, Campus Trindade, 

88040-900 Florianopolis, SC, Brazil 

Kafka, S.^ 

Department of Terrestrial Magnetism, Carnegie Institution of Washington, 5241 Broad 

Branch Road NW, Washington, DC 20015 

and 

Baptista, R 

Departamento de Fisica, Universidade Federal de Santa Catarina, Campus Trindade, 

88040-900 Florianopolis, SC, Brazil 

and 

Tappcrt, C. 

Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica, Vicuiia 

Mackenna 4860, 782-0436 Macul, Chile 

Received ; accepted 



^email: tiago@astro.ufsc.br 



-2 - 
ABSTRACT 

We report analysis of VRIJH photometry, and phase-resolved optical spec- 
troscopy of the eclipsing DA white dwarf plus dMe dwarf binary QS Vir. Model- 
ing of the photometric data yields an inclination oii — 74.9±0.6 and a mass ratio 
of g = M2/M1 = 0.50 ±0.05. Our Doppler maps indicate the presence of material 
in the Roche lobe of the white dwarf, at a location near the M star, likely due to 
accretion from the stellar wind of the M star (as opposed to Roche-lobe overflow 
accretion). We also constructed images of the brightness distribution of the M 
star at different epochs which reveal the location of two stable active regions. 
Doppler tomography shows that the majority of the Hydrogen and Ca II H&K 
emission originates on the active M dwarf, likely distributed in two preferred 
activity longitudes, similar to active regions on BY Dra and FK Comae systems. 

Subject headings: 
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1. Introduction 

The progenitors of close binary stars with a white dwarf as primary (i.e., more massive) 
component are thought to be initially wide detached binaries. In the course of its nuclear 
evolution, the more massive star expands first. If the separation is sufficiently small, this 
will yield a Common Envelope configuration, where both the remaining core of the more 
massive star and its lower mass companion revolve around each other within the expelled 
material from the primary. Responding to the enhanced angular momentum loss due to 
friction, the separation between the two stars decreases, and when the envelope material is 
expelled as a planetary nebula, the aforementioned close binary emerges. Such systems are 
known as post common envelope binaries (PCEB). 

Angular momentum is continuously lost via gravitational radiation and via a 
magnetically confined stellar wind from the main sequence star component {magnetic 
braking). As a result, the orbital period (and the binary separation) decreases with time, 
until the lower mass dwarf eventually fills its Roche lobe. Mass transfer then starts through 
the inner Lagrangian point (LI) of the binary, dumping material onto the white dwarf 
(WD) via an accretion disk or the magnetic field lines of the white dwarf. When the 
dwarf is a K/M star, those systems are known as cataclysmic variables (CV)Q The orbital 
period (and the binary separation) at which first contact occurs depends strongly on the 
radius/mass of the low mass main sequence star (secondary) and on the mass of the WD 
(primary). If the time required for the secondary star of a PCEB to come into contact with 
its Roche lobe is shorter than the Hubble time (~ lOGyr) the system is called a pre-CV 



(Schreiber & Gansicke 2003). 



^For a comprehensive description of CVs the reader is referred to Hellier (2001) and 



Warner (1995). 
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QS Vir (or EC 13471-1258) is an eclipsing WD+M dwarf binary system, discovered 



in the Edinburgh-Cape Blue Object Survey (Stobie et al. 1997). Kawka et al. (2002) and 



O'Donognue et al. (2003, hereafter DOD03) provided an extensive study of the system, 



combining photometric monitoring, time-resolved spectroscopy, and (in the latter case) 



UV spectra. By measuring the rotational velocity of the WD DOD03 shown that it is 



consistent with that of accreting CV primaries (Szkody & Mateo 1986). Therefore, they 



proposed that the system is not a pre-CV, but instead a hibernating CV; i.e. a former 



CV that became detached as the consequence of a nova eruption (Prialnik & Shara 1986 



Shara et al. 1986). Their photometric monitoring for over 10 years revealed small orbital 
period deviations from a linear ephemeris which they suggest is a result of random variation 
in the orbital period of the system, possibly caused by activity on the M secondary. A 
number of optical flares were observed, confirming that the M star is magnetically active. 
Phase-resolved optical spectra suggested that low level accretion is actually taking place in 
the system. 

We present here additional time-resolved blue spectra of the system using the RC 
Spectrograph at the CTIO 4.0m Blanco Telescope and Echelle spectra taken with the VLT, 
as well as recent near- infrared photometric data acquired with the Brazilian 1.6m Perking & 



Elmer telescope at Pico dos Dias. We also analyze the R and V(RI)c light curves of Kawka 



et al. 


(2002 


) and 


DOD03 



to explore the active regions on the M star, providing a pictorial guide to their location 
and evolution. This paper is organized as follows. In section [2] we present data acquisition 
and reduction procedures, while in section [3] we report the results from the analysis of the 
data with the different tools. We follow with a discussion of the combined analysis and a 
comparison with previous results in section |4j Final remarks conclude our paper in section 

El 
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2. Data acquisition and reduction 

2.1. Near-IR Photometry 

NIR photometry of QS Vir was collected with the CAMIV infrared camera attached 
to the 1.6 m Perkin & Elmer telescope of the Observatorio do Pico dos Dias (OPD/LNA), 
in southern Brazil, on 2008 March 18. CAMIV has a 1024 x 1024 pixels HAWAII detector 



cooled to 77K, designed for observations at 0.8 — 2.4yum (for further details see Rodrigues 



et al. 2006). The sky was patchy at the night of the observation and the humidity was 
very high, interrupting the observations before completion of a full orbit. Fast wind was 
also blowing at the dome, tilting the telescope in some cases. A log summarizing the 
observations is presented in Table [1} 



In order to avoid cycle to cycle brightness variations (very common in QS Vir, DOD03) 
the photometry was performed in a quasi-simultaneous mode, alternating between J and H 
filters after each dithering sequence. Data reduction was performed with CIRRED/IRAFQ 
The non-linearity, inherent to all reversed polarized semiconductor detectors, was corrected 



using the task IRLINCOR with the linear coefficient provided by Justiniano (2004). A 
mask for bad pixels was constructed using a normalized median of flat-field images. Sky 
level and dark frame subtraction as well as flat-fleld division and bad pixel mask corrections 
were applied to each frame. 

Aperture photometry was then extracted for the target star and 5 other nearby fleld 
stars, with the brighter fleld star being used as a reference star for differential photometry. 
This procedure minimizes the effects of sky transparency variations along the night and 



^IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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is particularly useful for non-photometric nights - such as that of our observations. In 
order to keep track of the changing position of the stars in the field due to the dithering 
procedure, we have used the task IMMATCH/CTIO. The relative magnitude of the target 
was converted to absolute flux (in mJy) using the absolute magnitude of the reference star 



together with zero point constants provided by the 2MASS project (Skrutskie et al. 2006). 



The data was phase-folded according to the mid-eclipse ephemeris of |DOD03 

Tmid(E) = HJD 2 448 689.64062 + 0.150757525 ■ E 
where E is the eclipse cycle. 



2.2. Spectroscopic Observations 

Most of our spectroscopic observations were obtained with the RC Spectrograph on 
the CTIO 4-m Blanco telescope on 2007 March 20. We obtained 46 spectra using the 
KPGL3 grating, yielding a spectral range from 3800-7500A at 1.2A/pixel, which provided 
a resolution of ~3A. A log of our observations is included in Table [!} For detector 
calibrations, we used the standard IRAF procedures, and for data reductions we used 
IRAF's ONEDSPEC/TWODSPEC packages. Due to non-photometric conditions (some 
passing clouds) we did not obtain spectroscopic standards, therefore our data are not flux 
calibrated. Example spectra are presented in Fig. [3j 

The Echelle spectrum was taken as a 30 min exposure on 2007 August 16 with UVES 
at UTl, ESO-Paranal. The DIC2 mode was used with cross-dispersers #2 and #4 and 
central wavelengths 3900 A and 7600 A. Reduction of the data was performed using the 
Gasgano interface with the ESO/UVES pipeline scripts. This included Bias subtraction, 
fiat-fielding, wavelength calibration, and the correction for the instrumental response 
function and atmospheric extinction using flux standard LTT 7987. The output consists 
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of three spectra with useful data in wavelength ranges 3800-4510 A, 5720-7480 A, and 
7700-9430 A, respectively. The spectra have slightly different resolving powers around 
70 000. 



3. Data analysis and discussion 



3.1. Photometry 



The V(RI)c data of DOD03 (see their Table 6) consist of four consecutive nights of 



observations in 1993. We will refer to the different sets of light curves from HJD 2 449 161 
to 2449164 as JDl to JD4, respectiveljf] and to our NIR light curves of HJD 2454 544 



(from 2008) as JD7. In addition, the R data from Kawka et al. (2002) consist of 2 cycles 
separated by approximately 3 days and will be referred to as JD5 and JD6. In Fig. [T] we 



present our JH light curves and the JD2 V(RI)c data set from DOD03 together with a 
model fit to the data. 

Fig. [l] shows that the eclipse of the WD (primary) becomes shallower and the 
ellipsoidal modulations of the secondary more pronounced with increasing wavelength, i.e., 
the contribution of the cool secondary becomes more important than that of the compact 
hotter primary at longer wavelengths. In order to probe the photometry for information on 
the system components we have applied a light curve synthesis code similar to the Wilson ^ 



Devinney (1971) code. For a detailed description of our procedure see Ribeiro et al. (2007). 



In short, we assume that the surface of the secondary star is given by one gravitational 
equipotential within its Roche lobe and we modify its radiation field to account for gravity- 
and limb-darkening effects. Because the primary is assumed to be a point-like source its 
ingress/egress features (partial eclipse phases) are not modeled. 



3JD1=HJD 2 449 161; JD2=HJD 2 449 162; JD3=HJD 2 449 163 and JD4=HJD 2 449 164 
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As a starting point to the fitting procedure, we estimated the WD fiux Fi at each 
wavelength from the depth of the echpses and assumed that Fi is constant throughout 
the orbit. We also obtained an initial estimate of the M star fiux F2 at each wavelength 
from the average fiux during eclipse. The best-fit solution is obtained by an algorithm that 
searches for the set of parameters which minimizes the normalized constraint function, 

where Fo{4>) and Fm{(p) are the observed and modeled fiuxes at phase 0, respectively, N is 
the number of points in the light curve and o"o(0) is the fiux error at phase 0. 



The light curves are fitted simultaneously for a single set of q, i and revalues. Since 
the relative contribution of each star changes with wavelength, we adjusted Fi and F2, the 
contributions of the WD and secondary respectively, accordingly in each data set. The 
resulting parameters (and corresponding values) are listed in Table [2] and were used to 
produce the model light curves of Fig. [TJ The values listed in Table [2] are the median values 
obtained through 10'^ Monte Carlo simulations, while the errors are the corresponding 
standard deviations. For each Monte Carlo simulation, the value of each data point in 
the light curves is changed according to a gaussian distribution with mean equals to the 
true data point and width equals to its photometric uncertainty. Each Monte Carlo light 
curve is then subjected to the same real data fitting procedure. Although a fraction of 
simulations (~ 20%) results in secondaries filling the Roche lobe, all resulting values are 
used to estimate the uncertainties. Table |2] list the radius of the secondary [7-2) and the 



volumetric radius of its Roche lobe (rRL, as proposed by Eggleton 1983), both in units of 
the orbital separation (a). We also show the resulting radius of the secondary in units of its 
Roche lobe radius (r2/rRL)- 



^The radius of the secondary star in units of the orbital separation (a). 
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The values obtained by the hght curve modehng procedure are in excellent agreement 
with those of DOD03[ which were independently obtained by combining the measurements 



of eclipse width and partial phases duration with Ki and K2 velocities measurements. 
However, the resulting 2) for the fit to the V(RI)c data is somewhat higher than 

expected for a good fit (~ 1.0), showing similar deviations in the (0-C) diagram for all 
filters (see bottom panels of Fig. [T]). The amplitude of the deviations in the JD2 light 
curves is similar in V and R, and lower by a factor of two in the I band. A similar analysis is 
performed for the JDl, JDS, JD4, JDS and JD6 data; the results are shown at the bottom 
panel of Fig. [!} Since the JDS and JD6 light curves are given in relative magnitudes (Kawka 



et al. 2002[ ), it is not possible to flux calibrate this data set. Therefore we use differential 



magnitudes for all our data in order to provide a direct comparison of the amplitude of the 
variations. Although different, the deviations from the model are present in all data sets. 
They are apparently larger for the JDl and JDS data and clearly larger for JDS and JD6. 

We checked whether the shape of the JD2 residual curve is a consequence of a wrong 
choice of system parameters. Given the symmetric shape of the ellipsoidal modulations, 
wrong values of orbital parameters {i, q and or stellar component flux {Fi and F2) 
would produce systematic residual peaks at quadrature (orbital parameters) or conjunction 
(component flux). The observed residuals are different from both predictions. Alternatively, 
a small error in binary phase (say, 6(f) ~ 0.02) could qualitatively explain the observed 
residuals. However, the amount of phase shift required to explain the residuals would result 
in a perceptible displacement of the sharp primary eclipses, which is also not observed. 
Therefore, the inferred systematic, phase-dependent residuals in the light curves seem to 
reflect intrinsic deviations from uniformity in the surface brightness distribution of the 
secondary star. We will return to this point later in the paper. We adopted the system 
parameters derived from the analysis of the JD2 plus JD7 light curves (which are flux 
calibrated) to compute the residuals of the JDS and JD6 data. 
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in Table [2] we list the flux for each component of the system in each filter. We 
used FITSPEC/SYNPHOT to search for the best-fit blackbody spectrum to the observed 
WD and LMD V(RI)c JH fluxes. We infer temperatures of Ti = 18000 ± SOOOiT and 
T2 = 2800 ± 200K for the primary and secondary, respectively. Furthermore, assuming 



radii of Ri = O.Olli?© and R2 = O.42i?0(|DODO3j) we find distances of rfi = 60 ± 20 pc and 
(i2 = 40 ± 5 pc. The large uncertainty obtained for the temperature and distance to the 
primary star is a consequence of using essentially red wavelength coverage to constrain the 
temperature of blue (hot) spectra. 



DOD03 discuss the presence of a flare in their JDS light curve at the beginning of their 



observations. Similar deviations from model light curves for a non-spotted secondary star 



were also reported by Kawka & Vennes (2005) (see their Fig. 1) with regard to the JDS 
and JD6 data. At the bottom panels of Fig. [l]we show an (0-C) of the JDS data where 
it is possible to check the resulting increase in flux around phase 0.25, when the system is 
close to quadrature. Variations of comparable amplitude are also seen in JD5 and JD6. 
Motivated by the amplitude and shape of these variations, we applied a light curve inversion 
technique ( [Lanza et al. 1998) to assess the location of the active region(s) responsible for 
the flaring (in the case of JDS) and the origin of surface brightness variations that could 
be responsible for the deviations of the JD2, JD5 and JD6 light curves from the model. 
Data from JDl and JD4 were not included in this analysis due to their incomplete orbital 
coverage. 



We have implemented the entropy inversion procedure described by Lanza et al. 



(1998), which has been successfully applied to data of active single M dwarf stars, in order 



to derive the position and number of spotted regions (Korhonen et al. 2007). Since the 



original description of the method concerns only single spherical stars, we included the 
tidal distortion of the surface of the M dwarf star of QS Vir. This allows our code to 



- 11 - 



reproduce the model light curves of FigjT] for a secondary star with uniform brightness. 
We then modified the intensities of a grid of surface elements, searching for the brightness 
distribution that yields the best fit to the data. The procedure attempts to simultaneously 
minimize the reduced of equation [2] and to maximize the entropy of the intensity grid 
(see Skilling 1987) in order to select the most featureless brightness distribution that fits 
the data. We must also keep in mind that the one-dimensional light curves do not provide 
information on the latitudinal position of spotted region and that, in this case, the entropy 
regularization algorithm tends to place structures at the center of the projected stellar disk. 

We show the resulting model light curves of the inversion process as solid lines in Fig. 
|2j Dashed lines shows the orbital modulation of an uniformly bright secondary star as a 
comparison. The gray scaled images presented in Fig. |2] are the changes in the surface 
brightness distribution of the secondary required to produce the observed deviations from 
the orbital modulation of a distorted star with uniform brightness. All maps are shown in 
the same gray scale to allow for a direct comparison of different filters at different epochs. 
Since broad band photometry measures essentially the photosphere of the star, positive 
deviations from the uniform brightness model light curve (dotted lines) represent active, 
hot regions and will appear bright in the images. On the other hand, negative deviations 
represent cooler regions and will appear dark in the images. 

The JD2 maps show that there are two different types of structures. Close to phases 
(f) ~ 0.2 and cf) ~ 0.7 we see two regions of increased emission that are ~ 0.5 in phase 
apart. It is also possible to see two regions of decreased emission at phases ~ 0.4 and 
~ 0.9, also ~ 0.5 in phase apart. These structures, which represent spotted (cool) regions, 
are present in all three light curves and are more evident at shorter wavelengths (V). An 
interesting aspect of the morphology of the light curve is seen when comparing the (0-C) 
for the data of JDl and JD4 to that of the JD2 and JDS data. Despite the incomplete 
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sample of the JDl and JD4 data, there is a strong resemblance between the amphtude of 
the residual of the JDl and JDS light curves and JD2 and JD4 light curves, suggesting that 
the brightness of the relevant active regions changes on timescales of day 

The results for JDS and JD6 are similar to those for JD2 in the overall shape and 
positioning of the features. We observe a small (A0 ~ 0.05) offset of the bright regions in 
the longitude of the two bright spots (from phase ~ 0.25 in 1998 to phase ~ 0.3 in 
2000). Furthermore, we notice that the brighter of the two spots is the one centered at 
phase ~ 0.8 in JD5 and JD6 whereas for JD2 the brighter spot is the one centered at 
phase ~ 0.25. Since the JD5 and JD6 data sets are not flux calibrated it is not possible to 
assess the correct flux level of the secondary star, which is essential for the fitting procedure. 
Therefore, caution must be exercised when comparing the results of the analysis of the 



DODOS and Kawka et al. (2002) data sets. Jin each individual data set it seems that active 
regions on this M star are stable and keep their location on time scales of, at least, a few 
days. "A comparison of the 1998 and 2000 data, suggests that these active regions may be 
relatively stable over much longer times, of the order of years. Two-spots configurations, 
0.5 phases apart are commonly present on single dM active stars, however their location on 



the surface of the star evolves with time (Vida & Olah 2006). What we witness in QS Vir 



is likely in accord with "permanent active longitudes" on RS CVn binaries and FK Comae 
stars, of active regions kept in place via tidal interactions of the two stellar components 



(Berdyugina & Tuominen 1998) 



3.2. Spectroscopy 

Fig. [sjpresents examples of our spectra, at four different orbital phases (inferior/superior 
conjunction and quadrature). The out-of-eclipse spectra are dominated by the underlying 
WD in the blue and the M star in the red. Ha appears purely in emission along with the 
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Ca II H&K lines, whereas the other Balmer lines show narrow emission cores on top of the 
broad absorption component from the WD. During eclipse, the only emission lines present 
are the Ha and the Can H&K lines, confirming that the latter and most of the former 
originates on the M star (Fig. |3l top). 



Fig. |4] presents trailed spectra for the Balmer and Can H&K emission lines and for the 
TiO molecular band. Apart from the Ha emission line, all other Balmer emission lines are 
single-peaked and follow the motion of the M star (Fig. |4]). Their K velocities are presented 
in Table [H 

The Ha line appears to have a more complex structure, presenting three distinct 
components. These components are also present in our Echelle spectrum of Fig. [5j In 
Table |4], we summarize the properties of the emission line components obtained after a 
decomposition of the line profiles into gaussians. The different component are shown as 
dashed lines in Fig. [5] and are labeled accordingly. This spectrum was taken at a different 
epoch, indicating that the observed line structures are either persistent or at least very 
common for the Ha line of QS Vir. The Echelle spectrum has a gap between 4500 — 5500A, 
therefore it is missing the H/? line. Notwithstanding, H7 shows the same configuration 
as Ha while only one extra feature is detectable on H(5. It is quite likely that the lower 
resolution of the CTIO spectra prevents us from distinguishing multiple components in the 
other Balmer lines (as indicated by the structure seen at H'^ in the Echelle spectrum). 
Possible mechanisms responsible for the asymmetric (and broadened) line profiles, is the 
Stark effect and flaring on the M star. Comparing our line profiles with the models of 



AUred et al. (2005) (also see Fig. 3 of Paulson et al. 2006) we favor a flare on the M star 



as the mechanism that causes asymmetric line profiles and variations on the EWs in the 
Balmer lines of QS Vir. 

The K velocity of the TiO band was measured by cross-correlating all spectra in the 
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data set with a template spectrum to find the best-fit velocity displacement in each case. 



As previously noticed by DODOS and Kawka & Vennes (2005), the spectral type of the 



template star affects the resulting RVs amplitudes and leads to systematic differences 
which are comparable to the internal errors of the fit. In order to avoid this, we have used 
the spectrum of QS Vir around mid-eclipse (0= 0.034) as our template. At these phases 
only the secondary star contributes to the observed spectrum as the WD is eclipsed. This 
method has the advantage of avoiding the uncertainties associated to the choices of spectral 
type for the secondary star, and determine the amount of line broadening one has to apply 
to a slowly-rotating template star. 

In addition to the selection of a correct template to correlate with the M star spectra, 
we have noticed that the choice of the spectral window has a strong impact in the 
resulting measured values. Specially, the selection of narrow spectral windows results in 
systematically lower values for the estimated K2, while broader ranges results in larger 
values. We proceed implementing an iterative procedure that searches for the spectral range 
that minimizes the residue of the template spectra and spectra of the star corrected for the 
resulting radial velocity. This procedure results in K2 = 256 ± 10 km/s with a wavelength 
range of A = [6940^5° : 7250l2o]A. The uncertainty in the estimative of K2 accounts for 
both the uncertainty in the fitting procedure and in the wavelength range selection. 



Our resulting K2(= 256 ±10 km/s) is consistent with the findings of DODOS, (K2 = 266 



± 6 km/s) and Kawka et al. (2002, K2 = 241±8 km/s). Adopting our inferred parameters 
(i, g, K2; see Table [2]) we find a binary solution with component masses of Mi = 0.6 ± O.S 
M0 and M2 = O.S ±0.1 M0. The computed binary parameters are listed in Table |2} 

For the two nights of our spectroscopic monitoring, we measured the EW of the full 
Ha line; the outcome is presented in Fig. |6] and the relevant measurements are presented 
in Table [5] The four spectra of night 1 were taken around phase 0.1 and have EW of the 
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full Ha line close to 6A; this is similar to the Kawka et al. (2002) quiescent Ha EW of 
5.6A. This measurement should represent a constant "basal" coronal emission from the 
M star, indicative of its activity level. On night 2, the EW of the line is highly variable, 
reaching ~-18A and -SA at phase 0.0 in two successive cycles. It is very likely that, at 
the second night of our observations, we observed a flare on the M star, with maximum 
strength at inferior conjunction of the star. The EW modulation is likely a combination of 
a projection effect (due to the rotation of the system) and flare evolution, decreasing its 
magnitude by almost 50% in the next cycle. Furthermore, the EW of the star declines to 
-6A between phases 0.5 and 0.75 but returns to a value of -SA at inferior conjunction. This 
abrupt "jump" of the EW value after phase 0.75 suggests that we either witness the end 
of the previous event or the beginning of a new one, best visible at phase 0.0. Taking the 
EW variation as an indication of the extent of the active region that is responsible for the 
observed activity, we assess that it should be centered at the back side of the M star. The 
EW measurements of the full Ha line are presented in Table |5| 



3.3. Doppler Imaging 



An alternative representation of the emission line structure can be provided using 
Doppler tomography (DT) techniques. It uses the changes in line profile with binary phase 
to compute a map of the line emitting sources in a two-dimensional Doppler velocity space. 
The technique assumes that (i) all points of the brightness distribution are equally visible 
at all phases; (ii) there are no intrinsic changes in brightness with binary phase; (iii) all 
movement is parallel to the orbital plane; (iv) all brightness sources are fixed in the rotating 



binary frame; and (v) the intrinsic width of the line profile is negligible (Marsh 2001). We 



apply the maximum entropy code provided by Spruit (1998) to produce Doppler maps for 



Ha, H/3, H7, B.6 and also for the Can H&K doublet. For a more detailed description of the 
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code, see Spruit (1998) 



In order to make a direct comparison of our results with DOD03 we have made the 
same assumptions with respect to the geometrical position of the WD and the M star. Since 
DT does not take eclipses into account, spectra taken during eclipse were excluded from the 
image reconstruction procedure. The results for all emission lines are presented in Fig. [7} 

The most prominent feature in all Doppler maps is emission component centered on 
the secondary star. For the Ha line, this feature corresponds to our component #2 of the 
trailed spectrogram (Fig. |4]). Our Ha map also presents emission from inside the Roche 
lobe of the WD (this corresponds to our component #1). The increased emission present at 
the back side of the secondary is the structure responsible for the component ^^3 or "Has" 
of the deconvolved spectra (see Fig. |4]and Table |3]). This indicates either that there is an 
additional source of emission at the back side of the M star, or that we observed the system 
when the M star was flaring, in agreement with our EW analysis conclusions of Fig. [6] 

A useful test to assess the validity of the results from DT is to compare maps produced 
for different half cycles (from phases to 0.5 and 0.5 to l.Oj^ this helps us to estimate 



what regions of the system are satisfying DT assumptions (e.g. Hoard et al. 1998 Marsh 
2001). The outcome of this exercise is presented in Fig. [s] The main difference observed for 
the first to second half cycle Ha maps, is a reduction in the intensity of the corresponding 
features #1 and #3 of the trailed spectra (Figs. Ill and [8|. The bottom panel of Fig. 



[8] shows the line profiles close to phases (j) = 0.25 and 0.75. If all assumptions of DT 
were satisfied, the two representations would be mirror images of each other (around zero 
velocity). However the low velocity component is almost absent in the second half cycle 



^Here phases and 1 are used only as reference. As we have stated earlier, spectra during 
eclipse were excluded from the imaging process. 
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{(f) = 0.75) implying that it is either fading along the orbit or that it has an anisotropic 
emission. 



4. Discussion 

It is challenging to pinpoint the source of origin and the mechanism for the various 
Ha emission line components. Component #1 is more intriguing, since it seems to change 
with time. Material inside the WD's Roche lobe was also present in the Doppler maps of 



DOD03 and was perceived as a persistent structure. Although the Doppler map of DODOS 
and ours are similai]^ the location of the emission component originating inside the Roche 
lobe of the WD appears to be different. Since we have only one orbit of the binary (thus 
we can not assess orbital stability) this feature is either fading during our observations 
(therefore it is transient) or it has an anisotropic irradiation field and is seen preferentially 



at phase 0.25 compared to phase 0.75. An inspection of the trailed spectra of DODOS 
(see their Fig. 15), reveals that the same structure is also much fainter close to ~ 0.75 
than ~ 0.25. However their Doppler map was constructed using all the data from their 
spectroscopic observations (thus data from different epochs), which likely has smeared out 
any information about possible time evolution of the feature and/or its location in the WD 
Roche lobe at different times. 

From our light curve models we derive a radius for the M star that is smaller than its 
Roche Lobe radius at the 2-a confidence level (see Table [2]). Therefore, accretion should not 
takes place through LI as a result of standard Roche lobe overflow of the M star. At the 
same time, wind accretion is a possible scenario to explain material that lingers inside the 
Roche lobe of the WD, likely trapped by its gravitational potential. Active chromospheres 



^with the exception of the lack of feature #S in the ^DODOS^ data 
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of single M stars reach ~1.1 stellar radii and are also easily blown up as a stellar wind 



(Donati et al. 2000). There is some evidence that mass loss via stellar wind of early-type 
M stars is of the order of 10~^^ 



]^g-i4 i\/r ,rv.-i 



Moyr- 



see 



Debes 2006 for a detailed discussion). 



The low temperature of the WD, compared to CVs of the same orbital period (Szkody & 



Mateo 1986), and lack of observational evidence for an accretion disc suggest that if there 



is accretion then it should be at a very low level. We conclude that a stellar wind should be 
the origin of the material inside the Roche lobe of the WD in QS Vir. Nevertheless, with 
such a low M it is hard to believe that any shock or collision will actually form from that 



material, which was the interpretation of DOD03 for the features inside the WD Roche 
lobe, also observed in our Ha trailed spectrum and Doppler maps. Material outside the 
Roche lobe or inside the lobe of the WD cannot be in a stable orbital motion with the 
system unless it is kept there by some factor other than gravity. An interesting alternative 
could be a gas cloud locked with the M star by a magnetic prominence system, as has been 



reported by Donati et al. (2000). 



Using a light curve inversion method we assess the details in the brightness distribution 



of the M star at six epochs of asymmetric light curves (using data from DODOS and Kawka 



et al. 2002). In Figure 2 we present our results for the epochs of observations that have 



good orbital coverage, namely JD2, JDS, JD5 and JD6. In JD2, there are two bright 
active regions centered at ~ 0.2 and ~ 0.7, which seem to be shifted slightly in phase 
(A0 ~ 0.1) in JDS. The JDS map, in turn, is dominated by a strong flare; which will be 
commented later. From the partial light curves of JDl and JD4, we have indications that 
these activity locations are also present during those epochs of observations. Interestingly, 
active regions at approximately the same region appears in JDS and JD6. 

The stability of active regions with time prompts to the consistency of active longitudes 



in RS CVn and FK Comae systems. Korhonen et al. (2009) find that such stable active 
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regions correspond to locations where magnetic flux tubes of opposite polarity emerge on 
the surface of the stars. Although similar high-resolution long-term data are not available 
for QS Vir, it will be interesting to investigate if such active regions evolve in latitude - 
as is suggested by the analysis of the JD5/JD6 data sets - in synchronism with activity 
cycles of the M star, or if such cycles have no effect on the location of the activity regions 
in QS Vir-type binaries. Possible implications on the level of magnetic braking in such 
systems (and consequently on the evolution of pre-CVs) remain to be determined. 

At JDS the flare appears during eclipse and disappears half a cycle later as the binary 
rotates. The flare is likely located between ~ 0.2 and ~ 0.4 (Fig. |2]), centered on 
~ 0.3. Furthermore, in our spectra, the EWs of the total Ha emission line indicate that a 
flare was also present during our spectroscopic observations. The EWs analysis place the 
source of the variation at the back side of the M star, which is conflrmed by the Doppler 
maps, reconstructed with the same data set. This does not come as a surprise, since this M 



star seems to be very active at all epochs of relevant observations (eg. DODOS). Regarding 
the decrement observed in the intensity of the Balmer lines - i.e. the component been 
stronger in Ha with decreasing intensity for the other Balmer components, we believe 



that the M star has experienced a solar-like flare ( Garcia- Alvarez et al. 2005). No similar 
structure at or close to that same region is found in the Doppler map of DODOS supporting 
the suggestion that this is a transient feature. 

A remarkable property of this system is its apparently high recurrence rate of flares 
which, in principle, could assist material flnding its way to the WD. In the context of 
close detached binaries, there is evidence that tidal interactions can induce the emergence 



of magnetic flux tubes at preferred longitudes on the surface of the star (Holzwarth & 



Schiissler 200S, for example). With this scenario it is possible to explain the presence of 



material co-rotating with the binary within the Roche lobe of the WD in QS Vir as material 
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trapped in magnetic loops from the M star. Changes of the position of the corresponding 
feature will occur accordingly with the regeneration/evolution of the gas cloud and the 



geometry of the magnetic field of the M star at that location. Donati et al. (2000) measured 
a flux ratio of Ha/H/? ~ 10 for the gas cloud of the post T Tauri star RX J1508. 6-4423; 
our respective value is Ha/H/? ~ 14 for the equivalent structure, in agreement with this 
interpretation. 



5. Summary and conclusions 

The main results of our time-resolved photometry and optical spectroscopy can be 
summarized as follows: 

• We provided light curve models of QS Vir, which (combined with the subsequent 
determination of K2) constraining of the binary parameters. The size of the M star 
yields a distance to the system of 40 ± 5 pc. We also find that the M star radius is 
smaller than its Roche lobe radius at a 2-a confidense level which, therefore, suggests 
that accretion does not take place via Roche lobe overflow in QS Vir. The main 
binary parameters results are summarized in Table |2} 

• Surface images of the M star were constructed and show four active regions with two 
different morphologies. There are two cool spotted regions 0.5 in phase apart and two 
hot active regions also 0.5 in phase apart. This configuration is persistent during a 
flare observed on the JD3 light curves and in the JD5/JD6 data. The different data 
sets JDl to JD4 and JD5/JD6 are separated by ~ 9 yrs. The spot responsible for the 
flare on the JDS data is centered at ^ 0.3 and occupies a sizable fraction of the star 
surface, with an azimuthal width of ~ 90°. 

• We have measured the RVs of the multi-component Ha emission line and of other 
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Balmer and Can H&K emission lines as well as of the TiO band. The results indicate 
that the chromospheric emission lines have larger K velocities than that expected for 
the centre of mass of the M star (as measured by the TiO band). The results are 
listed in Table [3l 

Analysis of our time resolved spectroscopy and Doppler images of the Ha emission 
line argue against the interpretation of the structures seen on the WD side of the 
binary as a result of low level accretion stream shock. It is likely that it is a result 
of wind accretion onto the WD. If the material inside the Roche lobe of the WD is 
long-lived, we suggest that it can be held in place by a prominence-like magnetic loop 



in a manner similar to that observed on single active stars (see Korhonen et al. 2009 
and references therein). 

The detection of Ha emission from the WD side of the binary was interpreted as a 
residual of a previous accretion state of QS Vir, suggestive that this is a hibernating CV 



(DODOS). With the data in hand, we can not examine the long term spectroscopic behavior 
of the system in order to compare the location and strength of the material inside the Roche 
lobe of the WD with time, and assess its stability and persistence. Our data also provide 
no indication that QS Vir bears characteristics of a hibernating CV. Since this is one of the 
few well-studied detached M dwarf- WD binaries, future observations complementing our 
work can assess the nature and stability of the active regions on the surface of the M star, 
and their role in angular momentum loss via magnetic braking, which drives the evolution 
of such systems. 
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Fig. 1. — The upper panel shows our NIR hght curves and JD2 data from O'Donognue et 



al. (2003) of QS Vir. Model light curves are shown as solid lines and different symbols are 



used for each filter. From top to botton, the NIR light curves are shown as open diamonds 
(H) and open circles (J) and the JD2 are in plus signs (I), crosses (R) and filled circles (V). 
The phases are repeated for better visualization. In the six bottom panels the residue shows 
deviations from the model (obtained for the NIR and JD2 light curves) to the VRI data of 
JDl, JD2, JDS and JD4 and the R data of JDS and JD6. The symbols are the same as in 
the top panel, for each filter. An horizontal line depicts level zero and vertical lines depicts 
phases of quadrature (0 = 0.25, 0.75,1,25 and 1.75). 
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Fig. 2. — Application of light curve inversion method on the differential light curves of 
JD2 (left), JDS (middle), JDS (upper right) and JD6 (lower right) data. The resulting light 
curves are plotted as solid lines, while the light curves for a uniform brightness distribution 
are shown as dashed lines. Below each light curve, we present the result of the light curve 
inversion method. Lighter regions are brighter. 
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Fig. 3. — Example spectra of QS Vir. At the top panel we show a spectrum taken during 
the eclipse, as can be noticed by the absence of the WD features . 
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Fig. 4. — Trailed spectra of the major Balmer lines (top four panels) and the Can H&K 
doublet and the TiO absorption band of QS Vir spectra. The position of the rest wavelengths 
for each line/band and radial motion of the secondary are depicted as solid lines. 
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Fig. 5. — Echelle spectra of Balmer and CanH&K lines of QS Vir are shown as solid lines. In 
order to increase the S/N we have combined the data in bins of O.IA. Dashed lines shows the 
resulting fit of Gaussians to the line profile with arrows indicating each component (according 
to the results in Table |4]). 
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Fig. 6. — EW of the Ha line (top) and Can H&K (bottom) vs time. At the top figure, we 
mark the orbital phases of the system. See text for details. 
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Fig. 7. — Doppler tomograms of Balmer and Can emission lines. Solid lines shows the Roche 
lobe for the parameters listed in Table [6] and crosses depicts the positions of the center of 
each component and of the center of mass of the system {v = 0). The black cross at the 
upper right panel of the Ha and CanK maps represents the corresponding velocity-space 
resolution of the data. 
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Fig. 8. — Upper panel: Doppler tomograms of Ha comparing the map produced using 
all data for the second night of observations (Pull) and for different half-cycles, of the same 
data set. A contour line depicts a region where the intensity is the same in all maps (set to 
the average of the middle panel map). Botton panel: Example Ha line profile separated 
by half a cycle. If all Doppler tomography assumptions were satisfied (see text for details) , 
the profile should be a mirrored image. The dotted vertical line depicts the zero-velocity 
position, i.e. assuming Ao(Hq;) = 6562. SA. 



Table 1. Log of observations 



Obs. type Date Filter No exposures Exp. time ndith Start End 



Photometry 18 03 2008 J 53 45s 4 2 454 544.591 2 454 544.659 

Photometry 18 03 2008 H 49 30s 4 2 454 544.595 2 454 544.660 

Spectroscopy 20 03 2007 46 180s - 2 454 177.833 2 454 177.846 

EchoUc 16 07 2007 1 1800s - 2 454 329.538 



Table 2. Model Parameters from the light curve fits and resulting binary parameters. 



The radius of the the secondary star (r2) and its volumetric radius {tlr^ Eggleton 1983) 

are given in units of orbital separation. 

Parameter V R I J H 

F2{mJy) 3.05 ±0.02 7.66 ± 0.04 28.2 ±0.1 68.3 ± 0.4 72.8 ± 0.4 
FiimJy) 3.42 ±0.02 2.60 ± 0.05 2.2 ± 0.2 0.9 ±0.1 1.0 ± 0.2 
2.74 2.62 1.20 1.05 0.94 

Simultaneously fitted parameters: 
i = (74.9 ±0.6)° g= (0.50 ±0.05) r^jj = (0.44 ± 0.03) 

r2 = (0.36+°:°^) r,/rRL = (0.82j:°f ) 

Binary parameters: 
Mi = 0.6±0.3Mq M2 = O.3±O.1M0 R2 = 0.43l^;E^^ R© 

d = 40±5pc a = 1.2±O.lR0 
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Table 3. Parameters of emission lines 



line 


7 (km/sec) 


K (km/sec) 


00 


Hai 


18.6±0.7 


158±18 


0.08±0.02 


Ra2 


-34.7±0.9 


307±8 


0.462±0.004 


Hag 


-23±2 


549±4 


0.457±0.004 


H/3 


-9.2±0.6 


307±5 


0.462±0.003 


H7 


-61.5±0.6 


297±5 


0.450±0.003 


m 


-20.9±0.9 


285±10 


0.467±0.006 


CaiiH 


53.4±0.2 


278±4 


0.452±0.003 


CaiiK 


14.2±0.2 


281±5 


0.451±0.003 


TiO 




256 ± 10 
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Table 4. Parameters of our Echelle spectra emission lines. 



Component 


V [km/s] 


EW 


Flux 


Ra - 1 


-158 


-5.19 


32.5 


Ha - 2 


51 


-1.25 


133 


Ha - 3 


226 


-3.00 


76.7 


H7 - 1 


-171 


-0.24 


2.86 


H7 - 2 


56 


-1.89 


21.9 


H7 - 3 


226 


-0.69 


6.97 


m- 1 


104 


-1.66 


19.0 


m- 2 


249 


-0.28 


3.30 


Can H - 1 


183 


-3.7 


49.4 


Can H - 2 


368 


-2.93 


6.90 


Can K - 1 


-116 


-0.25 


3.45 


Can K - 2 


119 


-1.75 


23.2 


Can K - 3 


253 


-0.69 


9.06 
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Table 5. EWs for the Ha line 



HJD phase EW [A] 



2454177. 


,8332 


0, 


.1039 


-5.75 


2454177. 


,8395 


0, 


.1452 


-5.82 


2454177. 


,8430 


0, 


.1689 


-6.30 


2454177. 


,8464 


0, 


.1913 


-6.15 


2454179. 


,7780 


0, 


.0042 


-20.44 


2454179. 


,7814 


0, 


.0267 


-17.01 


2454179. 


,7849 


0, 


.0494 


-17.64 


2454179. 


,7880 


0, 


.0700 


-18.04 


2454179. 


,7911 


0, 


.0906 


-17.67 


2454179. 


,7942 


0, 


.1112 


-18.11 


2454179. 


,7973 


0, 


.1318 


-17.11 


2454179. 


,8004 


0, 


.1524 


-17.55 


2454179. 


,8038 


0, 


.1753 


-16.26 


2454179. 


,8069 


0, 


.1959 


-16.04 


2454179. 


,8100 


0, 


.2165 


-14.76 


2454179. 


,8132 


0, 


.2371 


-14.16 


2454179. 


,8163 


0, 


.2577 


-14.31 


2454179. 


,8194 


0, 


.2783 


-13.65 


2454179. 


,8228 


0, 


.3008 


-14.77 


2454179. 


,8259 


0, 


.3214 


-16.46 


2454179. 


,8290 


0, 


.3420 


-15.67 
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Table 5 — Continued 



HJD phase EW [A] 



2454179. 


,8321 


0, 


.3626 


-14.74 


2454179. 


,8352 


0, 


.3832 


-14.69 


2454179. 


,8383 


0, 


.4038 


-13.11 


2454179. 


,8419 


0, 


.4279 


-11.94 


2454179. 


,8450 


0, 


.4485 


-11.40 


2454179. 


,8481 


0, 


.4691 


-11.21 


2454179. 


,8512 


0, 


.4897 


-9.71 


2454179. 


,8543 


0, 


.5103 


-9.33 


2454179. 


,8582 


0, 


.5361 


-9.39 


2454179. 


,8613 


0, 


.5567 


-8.49 


2454179. 


,8644 


0, 


.5773 


-8.42 


2454179. 


,8675 


0, 


.5979 


-8.52 


2454179. 


,8707 


0, 


.6185 


-8.31 


2454179. 


,8740 


0, 


.6409 


-7.89 


2454179. 


,8771 


0, 


.6615 


-7.33 


2454179. 


,8802 


0, 


.6821 


-6.21 


2454179. 


,8834 


0, 


.7027 


-6.41 


2454179. 


,8865 


0, 


.7233 


-5.51 


2454179. 


,8926 


0, 


.7640 


-5.89 


2454179. 


,8960 


0, 


.7865 


-5.50 


2454179. 


,8991 


0, 


.8071 


-6.04 
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Table 5 — Continued 



HJD phase EW [A] 



2454179.9022 


0.8277 


-6.34 


2454179.9053 


0.8483 


-6.38 


2454179.9084 


0.8689 


-6.74 


2454179.9115 


0.8895 


-6.88 


2454179.9196 


0.9432 


-6.66 


2454179.9232 


0.9672 


-10.28 


2454179.9267 


0.9901 


-10.07 


2454179.9301 


0.0132 


-8.68 


2454179.9336 


0.0363 


-8.14 



Table 6. Model Parameters used to calculate the position of the centre of mass of both 
components and to draw the Roche lobe on the Doppler maps. 



parameter 


value 


i 


75.5° 


MwD 


0.78 Mq 


M2 


0.43 Mq 


Porb 


0.15075 d 



